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According to Berry a wave-chaotic state may be viewed as 
a superposition of monochromatic plane waves with random 
phases and amplitudes. Here we consider the distribution of 
nodal points associated with this state. Using the property 
that both the real and imaginary parts of the wave function 
are random Gaussian fields we analyze the correlation func- 
tion and densities of the nodal points. Using two approaches 
(the Poisson and Bernoulli) we derive the distribution of near- 
est neighbor separations. Furthermore the distribution func- 
tions for nodal points with specific chirality are found. Com- 
parison is made with results from from numerical calculations 
for the Berry wave function. 

PACS: 02.50. Cw, 03.65.Bz, 05.45.Mt 

I. INTRODUCTION 

The nature of the quantum eigenstates in billiards, 
which are classically chaotic, has been subject to much 
theoretical and experimental work. The seminal studies 
by McDonald and Kaufmann pj of the morphology of the 
two-dimensional (2D) eigenstates in a closed Bunimovich 
stadium have revealed characteristic complex patterns of 
disordered, undirectional and non-crossing nodal lines. 
The spatial behavior of the eigenstates of chaotic billiards 
is still of considerable theoretical and experimental inter- 
est. For recent theory see, e.g., the review by Rob- 
nik f|, and references cited; examples of measurements 
on electron billiards and wave-dynamical analogues are 
found in, ejr (quantum corrals), &-[llJ] (microwave 
cavities), Jl^] (acoustic resonators), |13|] (surface water 
waves in tanks) and in |f4| in general. Other well known 
signatures of quantum chaos in closed billiards are re- 
lated to the distribution of nearest level separations and 
spectral rigidity. 

For open billiards, i.e., billiards with attached leads, 
the picture is less clear. One may use the poles of the 
scattering matrix which are related to the decay time 
from a billiard ]f5]dl7|]. When transport through a bil- 
liard takes place one may as an alternative focus on the 
fact that the wave function ip is scattering state with both 
real and imaginary parts. If we restrict ourselves to 2D 
systems, as we will do throughout this work, this means 
that there are two separate sets of nodal lines at which 
either Re[tp] or Im[ip] vanish. The intersections of the 
two sets at which Re[ip] = Im[ip] — define the nodal 
points. Numerical simulations have shown that the shape 



of distribution function for the nearest distances between 
these nodal points (DFNDNP) depends on whether the 
billiard is nominally either regular or irregular pj| . For 
transmission through chaotic billiards the DFNDNP ap- 
pears to have a general characteristic form, while for reg- 
ular billiards like, for example, rectangular ones there are 
specific features of the DFNDNP that depend on the par- 
ticular geometry, at least as long as only a few channels 
are open. Thus, besides the vivid physical role played 
by the nodal points as centers of vortical motion [^0] , 
their statistical distribution may tell if chaos is present 
or not. The present work relates to quantum transport 
transport in open electron billiards. The issue of wave 
function singularities is, however, part a much broader 
context ]2r^-|29[|. 

An appealing argument that favors our view that 
DFNDNP serves as a signature of quantum chaos is the 
coincidence with the corresponding distribution function 
for the Berry state pj| . According to Berry's conjec- 
ture a chaotic state can be viewed as a superposition of 
a large number of interfering monochromatic de Broglie 
waves |^(| 

i 

where Oj and 4>j are independent random variables and 
kj are randomly oriented wave vectors of equal length. 
The Berry wave function may be regarded as a standard 
measure of quantum chaos. In fact, there are beautiful 
experimental observations of Berry waves on the surface 
of water in an agitated ripple tank with stadium-shaped 
walls Q. 

So far all our conclusions about DFNDNP rest on nu- 
merical experiments JXSj ] . The Berry state is, however, 
available in a mathematical form that invites to analytic 
approaches. In the present work we will therefore model 
the DFNDNP f(r) and its main asymptotic behavior an- 
alytically using the fundamental property that the Berry 
function (1) is Gaussian random field pfjj. We will also 
show that there are other types of distribution functions 
that are related to the chirality of the nodal points. Each 
nodal point is a topological singularity of the wave func- 
tion ]2p| , pl| , p3| -p6t . As a result there is a vortex centered 
around each nodal point with definite chirality depend- 
ing on whether the current flows clockwise or anti-clock 
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wise as indicated in Fig, [|. 



II. DEFINITION OF VARIABLES 




Consider the Berry function ([!]) as the complex random 
function 



FIG. 1. Typical pattern of nodal lines Im[^(x,y)] = 
(black lines) and Re[ty(x, y)\ — (red lines) for the Berry 
function. Nodal lines in each set do not cross. Points at 
which the two sets intersect are the nodal points around 
which there is vortical flow in either clockwise (green dots) 
or anti-clockwise (blue dots) direction. 

We therefore label each nodal point by a = ±1. In 
analogy with /(r) we therefore introduce the distribu- 
tion functions f a ,a-'(r) for nearest neighbor separations 
between points with chiralities a and a'. Analytic ex- 
pressions for these distributions will be derived below 
and compared with numerical computations using the 
Berry wave function. As will be pointed out in the text 
our results partly overlap with recent work by Berry and 
Dennis ]26[ ] (the pair correlation functions g a ,<r' and the 
relation between the mean density and wave number k). 
Most recently Dennis has also considered the distribu- 
tion of nearest distances among nodal points using the 
Poisson model Q . 

In the following sections we will derive expressions for 
the distributions of nearest neighbor separations between 
points with chiralities a and a' . For this purpose we will 
also have to consider the pair correlation functions g a ,a'- 
We will make use of two different analytic approaches 
based on the Poisson and Bernoulli models. In addition 
we will also calculate the distributions by direct numer- 
ical methods, i.e., we locate the nodal points by simply 
computing the nodal lines for Im^) and Re(^f) and how 
they cross. In principle the numerical results represent 
the correct distributions and gives us a way to test the 
accuracy of the different analytic approaches. 



■i/;(r) = u(y) + iv(r) 



(2) 



where r is the 2D position with Cartesian coordinates 
xi,X2 and it(r), v(r) are two real random fields. We 
assume that w(r) and v(r) are mutually statistically in- 
dependent, homogeneous and isotropic Gaussian random 
fields with zero mean. The correlation function has the 
well known form 

a(s) = (u(r)u(r + s)) = (v(r)v(r + s)) = J (ks). (3) 

which is a direct consequence of the Berry function (Q) . 
To find the statistical properties of the nodal points Vj we 
have to consider the intersections of the zero level curves 
(nodal lines) of the fields u(r) and v(r), i.e., the roots of 
the two equations: 

ufo) = v(tj) = (rj-eR 2 ). 

As mentioned in the introduction the nodal points are 
the centers of vortices. The associated probability cur- 
rent J(r) is proportional to 

J(r) = Re [ip*(r) iV^(r)] = u(r)Vu(r) - u(r)Vti(r) . 

(4) 

More rigorously the nodal points are responsible for the 
vortices in the sense that the loop integral [pO 23 - 2j| 



<j> drv = <j> \78dr 

must be non-zero only if it encloses a nodal point. Here 9 
is the phase of the complex wave function. In the present 
work we consider the vorticity field 



u = VxJ. 



(•5) 



In our two-dimensional case it is normal to the (x, y)- 
plane, i.e., cj(r) = uj(r)h z where n 2 is the normal unit 
vector. Substituting (||) into (if) we have 

w(r) = [Vu(r) x Vu(r)] . (6) 

At the nodal point 

Uj = u(v 3 ) 

is the angular velocity of the current in the very vicinity 
of r j . In the following we will call uij the vorticity of the 
j-th nodal point. 
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III. GENERAL FORMULAS FOR THE DENSITY 
OF NODAL POINTS 



If we introduce the density of nodal points as 
K(r) = \ U { rj )\6(u(T))S{v(r)) 

we obtain 

3 

Let us introduce also another singular function 

K v (r)=w(r)5(u(r))S(v(r)) ^^(r-r,) 

i 

where 



(7) 



(8) 



(9) 



(10) 



equals ±1 for clockwise and anti-clockwise vorticities u>j, 
respectively. Therefore (|o|) defines the sign of the vor- 
ticity of the nodal points. Below we will refer to a as chi- 
rality. In ref. [^6[ it is named topological charge. There 
are as many points with a = 1 as with a = — 1. 
Next, let us define the mean density 

(K(v)) = p (11) 

and the correlation function for the random density 

G(s) = Wr)ft(r - s)) = {T S(r - r,)*(r - Tj - s)\ . 

\ ij I 

(12) 

Notice that because 7£(r) is a statistically homoge- 
neous and isotropic random field, the mean density p is 
constant and the correlation function G(s) depends only 
on the distance s between the points of observation r and 
r s. Formulas (Q), (g), ([□]) and @ form the basis of 
the statistical analysis of the nodal points distribution 
assuming that the functions u(r) and v(r) are random 
functions. 

The mean density defines a characteristic scale 



1 

7t 



(13) 



which we will use below as the natural unit of distances 
in the "gas" of randomly distributed points, i.e., we will 
use the dimensionless variable 



s 



VPS. 



(14) 



Moreover it is convenient to formulate some analytical 
results in terms of the dimensionless pair pair correlation 
function with dimensionless argument 



(15) 



We now introduce the mean density j(r) around of a 
given point. One can show that 
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7 (s) = - G(s) 



(16) 



which will play a crucial role in the following derivation 
of DFNDNP for the Berry function. 

We now consider some useful relations for the statistics 
of the nodal points. Firstly, consider the mean number of 
points inside a circle C r with radius r centered at some 
given point. It is obvious that the mean number of points 
enclosed by the circle is equal to 



(n(r)) = 2tt / j(y)ydy 



(17) 



Using the dimensionless coordinate (111) one obtains 



(n{£)) 



2tt I g(r)rdr. 
lo 



(18) 



This relation takes into account that the dimensionless 
correlation function g(r) is at the same time the dimen- 
sionless mean density. Secondly, consider also the rela- 
tion for the mean number of nodal points 



i(£)) = J2nP(n;e). 



Here P(n;£) is the probability that n neighboring points 
belong to the circle. These probabilities satisfy the nor- 
malization condition 



P(0;£) = l-^P(n;£) 



(19) 



n=l 



The probability P(0; £) is of great importance because 
it is directly related to the distribution function for near- 
est distances between a given point and its neighbouring 
points (DFNDNP) f m \ n {£). This is because the cumu- 
lative distribution of the nearest distances £ m in is given 

by 



PminW = P(4nn < £) = 1-P(0;i). 



(20) 



Therefore we may now write the following relation for 
the dimensionless distribution of nearest distances 



fn 



(21) 



Thus, the last formula reduces the problem of calculating 
the DFNDNP to that of finding P{0;£). Below we will 
find approximate expressions for P(0;£). However, we 
will first present asymptotic formulas for P(Q;£) and the 
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DFNDNP. For small I one may replace the exact relations 
(|l8|), ( ft9|) with the asymptotic forms 

(n(*))~P(l;*), P(0;*)~1-P(1;*) * - 0. 



Finally, from (|18|), (|2l| ) and the last relation above one 
obtains the folowing asymptotic formula for the exact 
DFNDNP 



fn 



2ir£g(£) [£ 0) 



(22) 



Let us now apply these general considerations to the 
Berry function (Q) . First of all we will calculate the mean 
density p (|ll|). Using the definition (Q) and the fact 
that the variables of the homogeneous Gaussian field and 
its derivatives are statistically independent at the same 
point we have 

p={\u>{v)\){8{u{v))){5{u{v))). (23) 

It is straightforward to show that 

1 k 2 

s , <M> = T 

where k is modulus of wave vector of the Berry function 

we obtain the 



The general behavior reminds superficially of the cor- 
relation functions for amorphous solids with short range 
order and distinct shell- like structures present. The cor- 
relations are, however, more long-range in the present 
case. The same pair correlation function was recently ob- 
tained by Berry and Dennis Pq] , although expressed in 
a different analytic form. The derivations are somewhat 
tedious as indicated by the Appendix. It is therefore re- 
warding that there is perfect numerical agreement with 
Berry and Dennis' results pr| . 

Using the expressions in the appendix we can find the 
asymptotic expression 



(6(u(r)))(6(v(v))) 



(24) 



fly). Therefore, substituting (J2J) into (|2_ 
final expression for the mean density 

(25 » 

This exact relation has been derived recently also by 
Berry and Dennis p6| 

Secondly, consider the density correlation function (|l2|) 



G{s) = (|w(r)w(r + s)\S(u(r))S(u(r - 



■s))*(«(r))J(«(r + s))) 
(26) 



The calculation of G(s) is given in the Appendix (All) 
and the dimensionless pair correlation function for arbi 
trary £ is plotted in Fig. 0. 





FIG. 2. The dimensionless correlation function g nodal for 
the Berry wave function ([!]) versus the dimensionless radius 
£. 



g{£) -1/4 (*->0) 



(27) 



from which we obtain the asymptote for the DFNDNP 
(p2|) at small £ 



(28) 



This exact result is useful for testing approximate ana- 
lytical solutions and numerical simulation data. 



IV. THE POISSON AND BERNOULLI 
APPROXIMATIONS FOR THE DFNDNP 

In order to model the DFNDNP by analytic means 
we may in a first attempt use the Poisson approxima- 
tion. This approach has been discussed recently also by 
Dennis p0| . The Poisson approximation implies that all 
■ points around a given one (which is located in center) 
are statistically independent, i.e., it neglects higher or- 
der correlations. Therefore we have to take into account 
correlations only between the given point and its neigh- 
bours. These correlations can be incorporated using the 
mean density of points 7 around the given point ( |l6| ) . Ac- 
cording to the Poisson law the probability that no other 
points belong to circle with dimensionless radius £ is 



P(0,£) = exp(-(n(£))) = cxp I -2tt J zg(z)dz 



(29) 



Using the relation (^lf), we easily obtain the formula for 
the DFNDNP in the Poissonian approximation 



/minM « 2%£g(£) exp I -2tt J zg(z)dzj . (30) 

One notes that for small £ the asymptote of the approx- 
imate DFNDNP @ coincides with exact one @. 

For the special case of uniformly distributed and com- 
pletely random points (g(£) — 1) we immediately obtain 
the well known result 133 JsT^l 



,/minW =27rtexp(-7tf 2 ) 



(31) 
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For convenience we also introduce the new dimension- 
less radius 



r = £/(£) 



(32) 



which refers to mean distance between nearest nodal 
points {£) . The main asymptotic for the DFNDNP (^|) 
then reads 



(33) 



A comparison between (|30|), (J33J) and the numerically 
calculated DFNDNP are given in Fig. [| 



A general disadvantage of the Poissonian approach is 
that all nodal points are competing with each other to 
be neighbours of a given point. It is clear, however, that 
only nearest neighbours of the given point actually par- 
ticipate in such a competition. Therefore we consider 
the Bernoulli approximation for the nearest distances of 
points which takes into account the competition between 
neighbouring points. Similar to the Poisson approxima- 
tion we again consider the circle C-ji of radius R with 
the center at a given point O and assume that all points 
except the given one are statistically independent. Fur- 
thermore, let us assume that the total number of points 
inside the circle C-n is just equal to the mean density 
integral 




FIG. 3. DFNDNP versus dimensionless distance r = £/{£) 
for (a) random points (|3l|) (dashed curve) and (b) for the Pois- 
sonian approximation (BO) (solid curve) with < £ >= 0.657. 
The straight dash-dot line is the corresponding asymptote 
(p3[). The histogram shows the distribution as obtained from 
direct numerical calculations of the positions of the nodal 
points Re('i/(xj,yj) — Ira(^{xj,yj) — for the Berry func- 
tion (|l|). In these simulations we have generated the nodal 
points in a large number of samples, typically of size (60 x 60) 
and with k = \/2n. The number of random plane waves in- 
cluded has ranged from 20 to 80. In the example shown we 
have included 40 plane waves and averaged over 200 samples. 
Except for statistical variations the same results are obtained 
also for other choices of the number of plane waves, sample 
size and value of k. 

Obviously the simplest model with (<?(£) = 1) cannot 
reproduce the true DFNDNP for the simple reason that 
the nodal points are not random points. The Poisson 
approximated function ( |30| ) is obviously in much better 
agreement with the numerical results although the dis- 
tribution falls off to quickly at large separations. The 
agreement for small z is more satisfactory with v = 0.765 
which is rather close to the value 0.68 obtained from the 
direct numerical calculations. Although the Poissonian 
modeling gives reasonable results we need to go beyond it 
for a better description of the intrinsic statistical, higher- 
order correlations among the nodal points as indicated by 
Fig. | 



n(R) = 2tt / g(£)£d£. 
Jo 



(34) 



With these conditions the distribution of each randomly 
located point belonging to C-jz point is exactly equal to 



n(R) 



(35) 



Next, let us consider another circle Ci with the center 
at the same origin O and radius I < R. Obviously the 
probability to find a point in this smaller circle is equal 
to 



p(£)= / f(£)d 2 £ 



c. 



n(R) 



where 



(n(£)} = 2ir g{r)r dr . 
Jo 



(36) 



(37) 



In the same way we have that the probability that a point 
does not fall into the circle Cg is equal to 



q(£) = l-p(£) = l- 



i(£)) 



n{R) 



(38) 



Since points are statistically independent the probability 
for all points to be outside the circle Ci equals 



P(0;£)= 1- 



x n[R) 



n{R) 



(39) 



From eqns. flgfj), (|21|), and (|3j) it follows directly that 
d f (n(£))^ n(R) ^ 1 

/m ^ ) = -P M = 2rfsm (i-LL| 

(40) 

To obtain the DFNDNP analytically we make the follow- 
ing approaches within the Bernoulli approximation, (i) 
In formula (^) we replace the number of points n(R) by 
the mean number of points (n(R)) (|37|). (ii) We choose 
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the radius R in such a way that inside the circle Cr there 
are a certain number of points which compete with each 
other to be the nearest neighbour to O. The minimum 
number of points in the circle is obviously three if we also 
include the central point. In the following we will simply 
use th is value. As a result we obtain, instead of formula 
( [AlCl) , the following expression 



/ min (*) £ 2ntg(£) ( 



1 



(41) 



It is easy to verify that this approximate distribution 
is normalized and has the same asymptote as the exact 
distribution function (|28|). 




FIG. 4. Plots of the DFNDNP for the Berry function (gj) 
versus dimensionless distance ( |32| ) with < I >= 0.658. Solid 
curve is the Bernoulli approximated distribution (p4|). The 
histogram is the same as in previous figure. 



In Fig. H the analytic results in (41) are compared 
with the numerical distribution obtained directly from 
the Berry function. The Bernoulli approach is evidently 
more powerful than our previous attempt in predicting 
the DFNDNP @. The reason would be that we now 
catch some of the higher-order correlations. 

The distribution ([ll]) has the same linear behaviour at 
small I as in previous expressions (|2|) and (|33|). The 
coefficient v = 0.678 which is quite close to numerics 
(0.68). 



V. THE MEAN CHIRAL DENSITIES AND 
CORRELATION FUNCTIONS 

To gain more detailed statistics of nodal points we con- 
sider statistical characteristics of chiral dependent nodal 
points density similar to (^) 



K, 



,(r)=^«5(r-r i J. 



(42) 



where j a numerates positions of vortices with chirality 
a = ±1. Formulas (g) and (0) can be written via the 
chiral densities (Ea) 



(43) 



Similar to (|Tl]) we introduce the chiral mean densities 

(K a (r)) = Pa = p/2. (44) 

The last equality follows from the mean spacial isotropy 
of the Berry function ([!]) which implies that in mean there 
is no preferred chirality of nodal points. 

Moreover we introduce the chiral the auxilliary corre- 
lation function 



(45) 



G v {s) = (K v (r)K v (v - s)) 

= (^2 a A r ~ r i)<rj5(r - r.j - s)^ . 

and correlation functions of the chiral mean densities (|4 
G^{s) = (7£ CT (r)ft CT ,(r-s)) 



From relations d 



= (Y,Hr-r l JS(r-v 3a ,-s)j. (46) 

dH), @ and © it follows that 

(47) 



G ua ,{s) = - [G(s) + aa'G v (s)}. 

Consider a nodal point with chirality a . Similar to 
we define the mean chiral densities around this point 



(It) 
(44) 



7<r<r'(s) = - Gw(s) 

p 



Then from (|47|) and (^8j) we obtain 



(48) 



(49) 



The subscripts (a, a') in ( [49|) indicate that the given 
point and its neighbours have the same vorticity, while 
(c, — a) refers to different chiralities for the given point 
and its neighbours. Therefore a knowledge of the cor- 
relation functions G(s) and G v (s) which are calculated 
in Appendix (formulas (A6) and (A24)) is enough to 



find the chiral mean densities and correlation func- 
tions. These chirality-dependent correlation functions 
are shown graphically in Fig. ^ in the scaled form 



(50) 



G 




FIG. 5. Plots of the dimensionless correlation functions 
Qa,a{i) = 5+1 9<r,-a(£) = ff-, and g(£) as denned in (jH^). 

Also in this case we find good agreement with Berry 
and Dennis |26|. Above we alluded to the general shape 
of pair-correlation functions for amorphous materials. 
With the dependence on a included we may obviously 
carry this naive picture a bit further to talk in a loose 
way about the nodal points as a two-component system 
in which 'objects' with the same vorticities (topological 
charges) repel each other. At the same time 'objects' 
with the same a may approach each other closely. Hence 
we may think about the system of nodal points as a thin 
slab of a fictitious binary amorphous solid or salt. Of 
course, this analogy should not be pushed too far. 

In order to find the distributions functions f aa >(t) for 
nearest distances between nodal points with chiralities 
a and cr' we insert g aa '{t) into (p0[). For the Bernoulli 
approximation one also has to label the mean number of 
points ([37]) as 



{n a ,<T>(£)) = 2-7T / g a ^(z)zdz 



(51) 



which is to be be inserted in (f40|) . We have found above, 
however, that already the Poisson approximation catches 
the gross features of the nearest neighbor distribution. 
For this reason and because the Bernoulli approximation 
is overly tedious we will restrict ourselves to the Poisson 
approximation at this stage. Thus we consider 



-2tt 



/o\<t(^) — 1 



l.oar 



(53) 



i 



The cr-dependent DFNDNP (||) for (a, a) and (cr, -rr 
are compared with numerical results in Figs |^ and 
The difference between the two combinations (cr, cr) and 
(cr, — cr) is very clear. In the first case there is a strong 
repulsion between the nearest neighbors, a result that is 
to be expected from Fig. |[ The distribution is, how- 
ever, of a very simple form. Essentially is corresponds 
a symmetric ring ('first shell') around the given point. 
Except for the low tail regions it is well approximated by 
a Gaussian, i.e., there is basically a random distribution 
within the 'shell'. 




FIG. 6. Distribution functions of nearest distances between 
nodal points with the same chirality a versus dimensionless 
distance (^) with < I >= 0.998. Solid line is calculated 
from function (|5^). The histogram refers to numerical work 
as described in Fig. |j| For this case we obtain the same value 
for < i >. 



(z)dz . (52) 




Using (|49|), ( |5ip and (|52j) one can show that the asymp- 
tote of f a ,-a(J) coincides with the asymptote for the 
DFNDNP (|28|). For j ' a , a {€), however, we obtain the quite 
different form 



FIG. 7. The same as in Figj^ but for nodal points with 
opposite chiralities a; < I >= 0.696. The histogram refers to 
numerical work as described in Fig. ^. For this case we find 
< I >= 0.726. 



7 



In the second case, (er, —a), the distribution shows how 
the nodal points are allowed to come arbitrary close to 
each other and how the shell structure is less pronounced. 
For small separations f a .-a is obviously the dominant 
term in the total distribution /. Finally we note that 
< I >~ 0.7 for opposite chiralities (topological charges) 
and < I >~ 1 for equal chiralities, i.e., there are 'in- 
ner' and 'outer shells' for opposite and equal vorticities, 
respectively. As shown by Figs |^ and [?] the Poisson ap- 
proximation reproduces the numerical results in at least 
a qualitatively correct way. 



VI. SUMMARY 

We have considered the distribution among nodal 
points associated with chaotic wave dynamics. The nodal 
points are of special interest as they are associated with 
vortical flow and chirality a which is cither +1 or — 1. In 
particular we have focused on the distribution of nearest 
separations among the nodal points. The reason is that 
distributions of this kind may carry information about 
chaotic dynamics as conjectured in |l9). We have intro- 
duced analytic approaches based on complex Gaussian 
random functions with the known correlation function 
Jo(fcs). Two cases have been considered, namely the 
Poisson and Bernoulli approximations. 

As a supplement to the analytic approaches we have 
performed numerical calculations to locate the nodal 
points and their vorticity using the Berry chaotic wave 
function in eq. ([j]). The numerical distributions com- 
puted in this way are in principle the correct ones and 
are therefore useful for testing the accuracy of the ana- 
lytic modeling. We thus find that already the Poisson 
approximation gives a good qualitative understanding of 
the distribution of nearest neighbor separations. Some 
higher order correlations are, however, not accounted for. 
For this reason we have also considered the Bernoulli ap- 
proximation which picks up some of these features. On 
the other hand the behaviour at large distances needs 
further work. 

When the distribution of nearest neighbor distances is 
separated into distances between nodal points with equal 
and opposite chiralities one obtains a picture that super- 
ficially reminds of binary amorphous matter. There are 
distinct 'inner' and 'outer shells' associated with oppo- 
site and equal chiralities, respectively. While points with 
opposite chiralities may get close to each other there 
is a strong repulsion among pairs with equal chirality. 
The reason is, loosely speaking, that nodal points with 
equal chirality have to appear in conjunction with 'anti- 
vortices' or saddle points. These 'anti- vortices' act as 
local 'beam splitters' and are necessary for two nearby 
vortices to spin the same way. These interesting fea- 
tures should be pursued further because the distributions 
among the different phase singularities are obviously in- 
terrelated. Effectively we may then arrive at a descrip- 



tion reminding of a ternary amorpous materials. 

The distributions discussed here relate to generic fea- 
tures of a wave-chaotic state. For this reason it would be 
of interest if they could be verified experimentally. Us- 
ing, e.g., micro- wave cavities this appears to be a real 
possibility pa]. 
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APPENDIX: CORRELATION FUNCTIONS FOR 
THE NODAL POINTS 

In this appendix we outline the somewhat tedious 
derivation of the various correlation functions. Consider 
the correlation function of the random density ( |l^ ) 

G(s) = 

(\u(t)u}(t + s)\5(u(r))S(u(r + s)6(v(r))6(v(r + s)) . (Al) 

Since u and v are statistically independent random fields 
it is sufficient to consider only the statistical properties 
of field u. First of all, note that the joint distribution 
of the values of the scalar Gaussian random field at the 
points r and r + s has the form 



W(u,u s ) = (S(u(r 
1 

exp 



u)8(u(r 



Us)) 



27iVl - a 2 (s) 



+ u 2 — 2a{s)uu s 



2[l-a 2 ( S )] 



(A2) 



Futhermore, we will need the reciprocal statistical prop- 
erties of Gaussian scalar field u and its gradient Vit. 
They are completely defined by the correlation vector 



e(s) = (u(r + s)Vu(r)) 



(A3) 



and the correlation tensor bij(s), i,j = 1,2 of vector field 
Vu. 

It is convenient to project this vector and tensor 
onto the local coordinate system related to the vector 
s through the longitudinal components 

"|| 0), u ||(r + s ) 
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and the transverse ones 

u±(t), u±(r + s). 

These components have the remarkable correlation prop- 
erties 

(u(r + s)u|| (r)) = e(s) , (u(r + s)ux(r)) = . (A4) 

The tensor correlation function becomes diagonal 

(u||(r)w||(r + s)) = 6||(s) (u±(t)u±(t + s)) = b±(s) , 

(A5) 



(u||(r)uj_(r + s)) = 0. 
Here we introduce the following notations 
d 2 a{s) 



{s)= _ bx{s) = _l^l, e(s) 

as z s as 



b ll (0) = b ± (0) = -k 2 = b. 



da(s) 
ds 



We may then write the density correlation function (Al) 

as 



G(«) 



1 



47r 2 [l-a 2 (s)] 



(\u\\v± - v\ l u±\\u\ ls v ±s - VpU± s \). 

(A6) 



For brevity we have here written the dependence on s as 
an index. The following averages are performed for the 
eight manifold Gaussian fields 

{«|| ( r )> W|| (r + s)}, {v\\ (r), v\\ (r + s)} 
{«j.(r),uj.(r + s)}, {v±(r),u ± (r + s)} (A7) 

Here each pair of variables has the correlation properties 



(u 2 ) = (v 2 )=bA(s), 

= K) = b, 

(u||W|| s ) = (wy«|| s ) = c(s), 
(UJ_U_L S ) = (v±v± s ) = bj_(s) 



(A8) 



where 



A( S ) = 



b[l-a 2 {s)]-e 2 (s) 



c(s) 



b[l-a 2 (s)} 
b ll (s)[l~a 2 (s)}~a(s)e 2 (s) 
b[l-a 2 (s)} 



(A9) 



It is convenient to transform to normalized random vari- 
ables 



Vbwy 



VbAisj' 



7T v± = 7=b 



With the correlation coefficients 

_ b ll (s)[l-a 2 (s)}-a(s)e 2 {s) 
a ~ b[l -a 2 (s)} -e 2 (s) ' 

expression ( |A6| ) now takes the form 

A(,s) 



(3 = 



bx(s) 



G(s) = p 2 



1 -a 2 {s) 



Ma, (3) 



where 



A(a,(3) = (\u\\vx - v\\ux\\u\\ B vxs - v^uxs 



(A10) 



(All) 



(A12) 



and the averaging is performed with respect to the dis- 
tributions 



w\\(u,u„) 
wx(u,u s ) = 



1 



27rVl — cr 
1 



: exp 



u 2 + u 2 



2cmu, 



: exp 



2(1 -a 2 ) 
u 2 + u 2 — 2(3uu s 
2(1 ~/? 2 ) 



Here distribution ion describes the statistics of "parallel 
pairs" {u||,U|| s } and {w||,«j| s } while wx describes prop- 
erties of "perpendicular" ones. 

It remains to calculate the function ( |A12| ) which will be 
done in two steps. Firstly we average over the statistical 
ensemble of parallel components {ui|,K|| s ,?;i|,tii| s } at all 
given perpendicular variables to obtain 



Ax(a,/?) = (|AA S |>_ 



(A13) 



where notation (•■•)_!_ means that all perpendicular vari- 
ables are considered as fixed, and that 



A = unvx - vnux, X s = un s vx s 



V\\ s UXs 



are two Gaussian variables with zero mean values with 
dispersions 



(Xl)=ul+vl 



and correlation 

(AA s )_l = a(uxUx s + vxvxs)- 
Let us normalize variables A and A s by the relations 
A A s 



±s 



J ±s 



which transform the function (A13) as 



Ax(a,0) = ^/(ul+vl)(ul s +vl s ){\zz s \)x (AH) 

where z and z s are Gaussian variables with unit disper- 
sion and the correlation coefficient as 



7 = {zz s )x 



V(«i + «x)(«L + «L)' 







Using the properties of Gaussian random variables one 
can derive after some algebra the following equation 



d 2 (\z lZ2 \)„ 
d 7 2 



A{5{ Zl )5{z 2 ))_ 
2 



Taking into account the initial conditions 



(A15) 



\Z\Z2\)X 



7=0 



2 

7T 



d{\ziz 2 \)_ 
(h, 



I 7=0 



we obtain from (A15) 



Hi) = (|«.|). 



2 r 



\/l — 7 2 + 7 arcsin 7 . (A16) 



Substituting ( |A16| ) into (|A14|) and averaging over the re- 
maining four perpendicular random variables we obtain 



A(a,/J) = (A ± (a,/3)) = (^K + + O^))- 

(A17) 

The angular brackets on the right hand of this equation 
imply an averaging over the _L-variables with the follow- 
ing joint distribution 



w(u,u s ,v,v s 



1 



4tt 2 (1-/? 2 ) 



x exp 



u 2 + v 2 + u 2 + v 2 — 20(uu s + vv s ) 



2(1 -P 2 ) 



1 



< I d£ drj £, 2 r] 2 exp 
'0 Jo 



1-/3 2 
( 2 + rj 2 - 2/3£r) cos \i 



2(1 -/3 2 ) 



(A20) 



If we use the new variables of integration p, d defined 
through 

£ = p cos (5 , r\ = p sin (5 



the integral (A20) transforms into the form 



2\2 



A( f i,P) = -(l-(3 2 ) 



d9 sin 2 ( 



dp p 5 exp 



ip 2 (l — /3cos/isin( 



/0 JO 

with 9 = 28. After integration over p we obtain 
A(ji,p) = (l-/3 2 ) 2 Q(/3cos/i). 



(A21) 



Q(z) = 



3z 



l + 2z 2 



(1-z 2 ) 2 (l-z 2 )5/ 2 



arctan 



(A22) 



Substituting (A21) into (A.19) we finally obtain 



(1-02)2 /-tt 

A(a,/3) = ^ — ---- / F(a cos n)Q((3 cos (j,)dfx (A23) 

27T ./__ 



In order to perform integration over four _L variables 
{u, u s , v, v s } we use the polar system of coordinates 

u± =£ cos 99, i?x=^sin^9 ii_L S = 77 cos -0, Uxs=??sin'^) 

which gives, instead of ( |A17| ). 

A(a,0) = (^(acos/j,)) (fi = ip-ip). (A18) 

Here the three random variables {£, 77, /1} are distributed 



where the function ^(7) is defined by th e equality ( A16| ) 
and the function (5(7) is given by ( A22 ). 

For the cr-dependent density correlation functions it is 
easy to show, using the Gaussian properties of the fields 
u and v, that the correlation function for cr-dependent 
random density @ has the form 

G v {s) = (R v (r)R v {r + s)} = p 2 - 2A ff - a{s)0(s). 



l- a 2(s) 



(A24) 



2tt(1-/3 2 ) 



exp 



£ 2 + rj 2 — 20^-qcosii 



2(1 -P 2 ) 



£, r) > 0, /US [-7r,7r]. 

To exclude the random variables £ and 77 we rewrite re- 
lation (A18) as 



A(a,/3) = ±- C T{ 
2tt J 



a cos n)A(fi, 8)d[i 



(A19) 



where 
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